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ABSTRACT 



minor axis. By using the X-ray data in tandem with available radio and opti- 



r \ • Although rapid reorientation of a black hole spin axis (lasting less than a 

_H ' few Myr) has been suggested as a mechanism for the formation of wings in X- 

Ph' shaped radio galaxies (XRGs), to date no convincing case of reorientation has 

O ■ been found in any XRG. Alternative wing formation models such as the hydrody- 

-^ ' namic backflow models are supported by observed trends indicating that XRGs 

, ^, ' form preferentially with jets aligned along the major axis of the surrounding 

medium and wings along the minor axis. In this Letter, we present a deep Chan- 

^ ■ dra observation of 4C +00.58, an oddball XRG with its jet oriented along the 

OO 
(N 
G^ ■ cal data, we estimate relevant timescales with which to evaluate wing formation 

. ■ models. The hydrodynamic models have difficulty explaining the long wings, 

(^ ■ whereas the presence of X-ray cavities (suggesting jet activity along a prior axis) 

O . and a potential stellar shell (indicating a recent merger) favor a merger-induced 

reorientation model. 

>< : Subject headings: gal™: active - gala^es: individual (4C +00.58) 

1. Introduction 



X-shaped radio galaxies (XRGs) are double-lobed radio galaxies ( iLeahy fc Williams 



19841 ) which also possess a pair of long, faint, centro-symmetric "wings". They have gained 
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notoriety as a possible signature of a rapid (within a few Myr) reorientation of the supermas- 
sive black hole (SMBH) spin axis, conceivably induced by galaxy mergers in which either 
accre tion torque or a SMBH merger causes a spin-flip (JRottmanru l200ll : iMerritt &: Ekers 
20021 ). In this scenario, wings are "fossils" tracing the prior jet axis which radiatively decay. 



However, no convincing case for reorientation has been made in any individual XRG, 
whereas several lines of evidence support a hydrodynamic origin for the wings. For instance, 
in most XRGs, the wings are coaligned with the minor axis of the host galaxy and the jet 
with the major axis ( Capetti et al.ll2002l : ISaripalli fc Subrahmanvan 20091). A similar major- 
axis — radio alignment trend exists in the X-rays ( iHodges-Kluck et al.ll2010l ). These results 
have been ii iterpreted to favor models in which the wings are produced by r adio- lobe — gas 
interaction dheahv fc WiUiamslllQSJ : IWorrall et al.lllQQsi ICapetti et al.ll2002l ). Additionally, 
no clear signs of mergers h ave been found in a spectroscopic study of XRG hosts, whereas 
they may be overpressured (JLandt et al.ll2010l ). 



In this Letter, we identify 4C +00.58 (Fig. 1, classified as a candidate XRG by I Cheung 



20071 ) as one of the best candidates for a merger-induced reorientation based on quantities 
derived from a deep Chandra X-ray Observatory observation and publicly available data. 
Unlike other XRGs, the 4C +00.58 jet is coaligned with the minor axis of its host. Even if 
only a fraction of XRGs are produced by reorientation, their frequency may be important 
for estimating Laser Interferometer Space Antenna detection rates. 



We use a Galactic absorption of N^ = 7.14 x 10 cm (iKalberla et al.ll2005[ ). as well as 



the Wilkinson Microwave Anisotrop y Prohe cosmology {Hq = 71 km s ^ Mpc ^, ^a 
and Qm = 0.27; ISpergel et aDl2007h . At a redshift of ^ = 0.059, 1" = 1.13 kpc. 



0.73, 



Observations 



We obtained a 93 ks Chandra exposure toward 4C +00.58 using the Advanced CCD 
Imaging Spectrometeiij (ACIS) an d combined it with a 10 ks archival observation (obs. IDs 
10304 and 9274; latter published in lHodges-Kluck et al.ll2010l ). The source is centered at the 
nominal aim point on the ACIS-S3 chip. The data were reduced with the Chandra Interactive 
Analysis of Observatio ns (CIAO v4.0) software, and spectral fitting was performed with 
XSPEC (lArnaudlll996l ). We extracted a 0.3 — 10 keV bandpass lightcurve (binned to 600 s) 
from empty regions to check for background flares, but found no 3a deviations. The extended 
emission around 4C +00.58 is less than 45" in radius, so we use local background for spectral 
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extraction. 



We also use a NRACp Very Large Array (VLA; iThompson et al.lll980l ) 1.4 GHz map 
( Hodges-Kluck et al. 120101) and a 4.9 GHz map produced by combining archival snapshot A- 
array data from lBest et al.l (119991) and C-array data from program AC406. We use Sloan Dig- 
ital Sky Survey (SDSS: lAdelman-McCarthy et al.ll2008l ) red (623.1 nm) and green (477.0 nm) 
images from two 54 s exposures of the host (SDSS J160612. 68+000027.1) to measure color 
and magnitude, correcting for the smaller point spread function of the green images as well 
as sky background, the 1000-count software bias, and Galactic extinction. 

The 1.4 GHz map is shown in the left panel of Figure 1. The primary lobes of the radio 
galaxy lie nearly on an east-west axis and have a well defined boundary, whereas the faint 
wings are oriented in a north-south direction. The jet experiences a dramatic bend (by 60°) 
just before terminating, and the cocoon in the 1.4 GHz map is notable for a well defined 
edge with a surface brightness about five times that of the wings. The 5 GHz map resolves 
the jet into a string of knots (§3.3). We detect no counterjet. 

The X-ray emission (Fig. 1) is made up of two components: bright emission spatially 
associated with the jet/AGN and a compact diffuse atmosphere. To isolate the atmosphere, 
we mask point-like sources and restrict the energy bandpass to 0.3 — 3 keV (Fig. 2); this 
energy range contains 80% of the photons within 45". On the basi s of an unsharp mask 



image and a weighted Vq ronoi tessellation adapatively binned image (JDiehl &: Statlerl 12006 



Cappellari &: Copinll2003l ). we have identified several X-ray cavities (Fig. 2). These cavities, 
labeled Cne, Cnw, Cse, and Csw, are low surface brightness regions in the hot atmosphere 
bounded by "spurs" of greater surface brightness. The cavities are deep negatives in an 
unsharp mask image in which a heavily smoothed map (40 px) is subtracted from a lightly 
smoothed (3 px) map. The residuals are shown in Fig. 2. 

The SDSS image (Fig. 3) reveals a dim, resolved extension to the southwest of the host 
galaxy. The elliptical host shows no obvious internal structure, but th e extension m ay be 
a stellar shell from a prior minor merger (for a discussion of shells see iQuinnI Il984j ) . The 
extension is red {g — r ~ 0.7), but bluer than the host {g — r ~ 1.0). The apparent magnitude 
of the entire extension is m^ ~ 20. 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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3. Results 

3.1. Hot Atmosphere 

The bright X-ray emission immediately near the jet is nonthermal, otherwise the ex- 
tended emission is consistent with isothermal plasma out to 50 kpc from the AGN, with 
kT = 1.1 ± 0.2 keV within 25 kpc and kT = 0.9 ±0.1 keV outside. The photon statistics 
preclude a deprojection analysis. On the basis of the resulting emission measures, densities 
within 25 kpc are a fewxlO^'^ cm~'^, implying pressures P = 10~^^ to 10~^^ dyne cm~^. 
Thus, we adopt a sound speed Cg ~ 400 km s~^ throughout the region. The total luminosity 
of the hot atmosphere is Lx ~ 3 x 10^^ erg s~^. 

The gross morphology of the X-ray e mission does not coincide with the host galaxy. 



Using only the 10 ks exposure, we argued ( JHodges-Kluck et al.ll2010l ) that the orientation 



agreed with the host although the elhpticity did not. With deeper data it is evident that 
while we correctly excluded nonthermal emission to the southwest, we could not exclude 
nonthermal emission to the northeast. 

The most notable feature of the diffuse X-ray map is Cnw (Fig- 2), which is coUinear 
with CsE and the AGN. Cnw is enclosed by emission and cospatial with a spur in the radio 
cocoon, suggesting that Cnw and Cse are jet-blown cavities. Cne and Csw are also collinear 
with the AGN and are associated with the bases of the wings, with walls extending into the 
surrounding medium. 

The free-free cooling time of the Cnw and Cse bounding material places an upper limit 
on their ages: 

^ff ~ It%^ ^ 1-6 X ^0%'^'nZl JT ~ 500 Myr (1) 

where T7 ~ 1.4 is the temperature in units of 10^ K and n_3 ~ 4 is the density in units of 
10~^ cm^^. If Cnw is a spherical bubble (of radius 5 kpc) inflated at a locally estimated 
pressure P = 5 x 10~^^ dyne cm~^, the work done to inflate it is PF ~ 7 x 10^^ erg. Since 
the transonic expansion time is 10 Myr, the minimum average kinetic luminosity of the jet 
during that period is 2 x 10"^^ erg s~^ if Cnw is jet-blown. 



3.2. Wings 

The long wings are associated with low-signal X-ray structure (Fig. 2, bottom right). 
The approximate wing symmetry implies coherent formation, but the southern wing is 
slightly longer. 
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If the wings expanded transonically, their length (~ 36 kpc) imphes an age of 90 Myr, 
although this is a minimum age since the wings are seen in projection and may expand 
subsonically. Conversely, the synchrotron decay time tgync provides a maximum wing lifetime 
assuming the radio emission traces the entire wing v olume and the wings were inflated 



by the radiating plasma. To estimate tgync, we follow iTavecchio et al.l (120061 ) and use the 
1.4 GHz map to estimate the equipartition field B^q. We take the spectral index, a = 0.7 
(§1, oc z^~°), from low-resolution radio photometry and assume 7min ~ 10. For cylindrical 
wings of r = 6 kpc and h = 36 kpc, we obtain B^q ~ 10 fiG. We then find the electron 
Lorentz factor 7 ~ 5600 from UsiiA ghz] = 4 x 10^^i?7^ = 1.4 GHz, and thus find 

isync ~ 2.4 X 10^4-^5;^^ yr ~ 40 Myr (2) 

where 74 is in units of 10^ and B^q is in fiG. This value represents the cooling time of the 
1.4 GHz electrons; we emphasize our assumption that this is "first-generation" plasma occu- 
pying the wings, tsync is insensitive to projection effects relative to the transonic expansion 
time: if the wings are longer by a factor of 2, B^q decreases by a factor 2^/*^'^'*'") ~ 1.2 and 
tsync increases by a factor ~ 1.4. The disagreement between tgync and the expansion time 
suggests either supersonic expansion (i.e. like a jet-blown cocoon) or wing replenishment by 
supersonic inflowing lobe plasma. 

The work required to inflate the wings, assuming the cylinders above, is PdV ~ 10^^ erg. 
An age of 40 Myr implies an average kinetic luminosity of Lkm ~ 8 x 10'^^ erg s~^ applied to 
the wings alone. 



3.3. Nucleus and Jet 

The strong central X-ray point source corresponds to the AGN and contains 1100 counts. 
The spectrum is flt well by a model consisting of an unabsorbed power law with spectral 
index a = 0.7 ± 0.1 and a weak thermal model frozen at kT = 1.0 keV. The nonthermal 
luminosity is Lx ~ 7 x 10^^ erg s~^ between 0.3 — 10 keV. 

The X-ray jet is the next brightest feature and traces the radio jet well, including several 
X-ray knots (Fig. 4). A super-sampled X-ray image reveals that two of these knots coincide 
with more compact knots visible in the 5 GHz image and a bright region in the 1.4 GHz 
jet. There is a modest radio (1.4 — 5 GHz) spectral gradient between the inner and outer jet 
(Table 1), and the outer jet does not appear in the X-rays. Along the inner jet, we measure 
a broad-band radio to X-ray spectral index, arx ~ 1-0, which is consistent with the X-ray 
spectrum (a^ ~ 1.1 ± 0.4). Thus, the X-ray emission is consistent with a synchrotron origin 
(requiring a concave-down spectral energy distribution). Alternatively, assuming an inverse 
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Compton origin for the X-rays, we follow iTavecchio et al.l ( 120061 ) to estimate the jet Doppler 
factor 6 = [7(1 — /3cos^)]~^ in the knots by finding B such that -Beq^ = Biq/5. This results 



in (5 ~ 10 (Table 1) and -Bknot ~ 10 /iG. As the jet X-ray emission is unlikely dominated by 
inverse-Compton emission, 5 is formally an upper limit. -Bgq is derived with no knowledge 
of the X-ray emission and predicts P ~ 3 x 10~^^ dyne cm~^, in agreement with spectral 
fitting estimate. 

Although no counterjet is visible, the 1.4 GHz emission suggests the jet has point sym- 
metry. If the jet tail is dragged, the tail will cool radiatively with its length set by the cooling 
time (i.e. no additional plasma influx). We derive an upper -Beqi^ ~ 120 /xG at z/^ = 5 GHz 
and thus find tgync ~ 0.5 — 15 Myr for 5 = 1 — 10. We assume 5 declines along the jet, so for 
(5^5, tsync ~ 6 Myr. 

The projected length of the radio cocoon gives a transonic expansion time of 35 Myr. 
The work done to excavate a cocoon with semimajor axis a = 17 kpc and semiminor axes 
6 = c = 8 kpc is PdV ~ 7 x 10^^ erg, so the minimum average kinetic luminosity while 
inflating the cocoon is L\an ~ 7 x 10^° erg s~^. 



4. Wing Formation Models 

We consider three wing formation scenarios: an overpressured outburst, conical preces- 
sion of the jet axis, and merger-induced reorientation of the jet axis. 

In the backflow models, wings are produced by pressure- or buoyancy-driven back- 
flowing plasma from the terminal shocks evolving in the hot medium. Since we see no 
obvious evidence for a plumed jet directly feeding the wings, the most pla usible of these 



mode ls is a "blow-out" from an overpressured cocoon early in the source's life ( ICapetti et al. 



2OO2I ). In this model, the native atmosphere confining the young jets is aspherical with a 



preferential direction along which the cocoon ruptures. 

The backflow model suffers from several difficulties. First, the wings must expand at 
most transonically. The long projected length of the wings in 4C +00.58 (requiring an 
AGN lifetime of at least 90 Myr) is difficult to reconcile with the cocoon length unless the 
cocoon inclination angle from the line of sight 6^1, ns is less tha n 30°. The equivalent widths 



of O II (37 27A) and O III (5 007A) measured by iLandt et al.l ( I2OIOI ) argue against a steep 



inclination ( Landt et al.ll2004j ). Furthermore, the maximum wing lifetime tsync disagrees with 



the transonic expansion time, although the wing plasma may be continu ously replenished. 



Secon d, the model relies on strong backflows typically associated with iFanaroff fc Riley 



fll974j . FR) type II sources. 4C +00.58 is not easily classified as FR I or II, but at Mr 22.7 
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and logLr(1.4 GHz) ~ 25.3 W Hz ^ falls very clo se to thelLedlow fc OwenI (Il996l ) boundary 



(-^radio oc -^opt). Smce XRGs generally l ie near thelLedlow &: Owed (11996) FR I/H boundary, 



they may be a transition population (ICheung et al. 



20091: 



Landt et al. 



20101 ). Finally, the 



radio galaxy is misaligned with its host, so the ICapetti et al.l ( 120021 ) model cannot produce 
wings until the jet escapes the interstellar medium. 

Conical precession is a simple model in which the jet axis swings around, so the wing 
extensions and post-bend jet are equidistant from the AGN. This requires a steep inclination 
angle of 6'los < 30°. In this model, the wings trace the jet history, and the oldest plasma 
can be no older than the synchrotron cooling time. To obtain the precession rate (j), we 
compare the projected length of the post-bend jet to its cooling time (about 6 Myr) and 
obtain ~ 4° Myr~^ (corresponding to a mildly supersonic v ~ 460 km s~^). A 180° rotation 
takes ~ 45 Myr and implies supersonic expansion for the cocoon. The timescale is consistent 
with the wing tsync, but the model must also explain the larger far-side cone, as both the 
far-side cocoon and wing have longer projected lengths. Nota bly, precession does not explain 
the presence of Cnw and Cse, and numerical simulations (e.g. iFalceta-Goncalves et al.ll2010l ) 
indicate that it would not preserve obvious cocoon structure. 

In the reorientation (spin-flip) scenario, the wings are fossil lobes of a jet whose direction 
rapidly changed, either due to accretion torque or coalescence of a SMBH binary. There is 
circumstantial evidence for such a spin-flip: a possible stellar shell indicating a minor merger, 
and the cavities Cnw and Cse with overlapping cocoon extensions implying somewhat recent 
jet-gas interaction along an old axis. We describe a reorientation scenario for 4C +00.58 
presently. 

Given the small size of Cnw and Cse, the jet was in a weak or "off" state prior to the 
minor merger, but the SMBH spin axis was aligned with the major axis of the host. Upon 
ignition, the jet quickly formed Cnw and Cse- However, since the angular momentum axis 
of the accreting gas is generally misaligned wit h that of the SM BH, accretion torque will 
reorient the black hole's spin within a few Myr ( iDotti et al.l 120091 ). Hence, Cnw must have 
expanded at Vexp > 2cs if no prior cavity existed. 

Once accretion torque moved the jet to the wing axis, it inflated the wings as active 



exp 



< 



lobes. The transonic lateral expansion time of the wings is 10 Myr, implying v, 

during wing inflation. The "Z-shaped" wing extensi ons (covering 30° ) could be explained 



either as post-reorientation "wiggles" from a hot disk ( jDotti et al. 



2009, show that a hot disk 



is required for wiggles o f this magnitude) or intera c tion betwe en the lobes and merging ISM 
swirling into the host ( JGopal-Krishna et al.l l2003l : IZierll2005l ). These Z-shaped extensions 
then ev olve buoyantly, and may b e replenished by the primary lobes once the jet axis has 
moved ( IGopal-Krishna et al.ll2003l ). We suppose the jet moved to its current position due 



to coalescence of the SMBH binary or ongoing accretion torque, then formed the present 
cocoon. Since the jet may have experienced small realignments during the wing inflation 
phase, we infer a reorientation timescale of fewer than 50 Myr since jet ignition, well within 
the free-free cooling time of the Cnw walls and constrained by the wing tsync- Although it 
is possible that the system represents only a single spin-flip from the wings to the present 
location, this hypothesis does not not explain Cnw or Cse- 

The presence of a stellar shell must be confirmed, and it is possible that Cnw and 
Cse are not jet-blown cavities but rather part of a cocoon-evacuated shell (with bounding 
material describing a ring perpendicular to the jet) produced by a jet-ignition shockwave. 
The overlap of Cnw by the radio cocoon is then due to backfiow filling the cavity. Assuming 
a circular ring, the eccentricity of the ring implies 6'los ~ 60°, far above the 30° required to 
reconcile the cocoon and wing lengths. 



5. Context and Summary 

There are few deep X-ray observations of XRGs. Apart from 4C +00.58, there is a 



~ 100 ks C/ianc?ra observation of NGC 326 and a 50 ks exposure of 3C 403 (iKraft et al.ll2005l ). 
These data show X-ray emission on different scales and of differing surface brightness, forcing 
the hydrodynamic hypothesis to contend with a variety of environments. The cavities in 
4C +00.58 also demonstrat e that X-ray observations a re not only useful for studying backfiow 



models. Our prior survey (JHodges-Kluck et al.l 120101 ) and this study suggest exposure times 



of at least 100 ks are required to examine detailed structure. 

We know of no clean evidence for merger-induced reorientation. Even in our toy model, 
the wings are produced by merger-induced accretion rather than a n instantaneous spin-flip , 



so the black hole merger itself would involve mostly- aligned spins ( iBogdanovic et al.l 120071 ). 
Nonetheless, the presence of an apparent stellar shell suggests that searching for structure 
in the hosts of XRGs may provide strong indirect evidence for such mergers in a subclass of 
these objects. 

We have presented a deep Chandra observation of the XRG 4C +00.58. The hot 
atmosphere is roughly co-spatial with the radio galaxy and has a temperature of kT ^ 
1.0+0.2 keV. An X-ray jet of about 5 kpc is detected, overlapping well with the 5 GHz knots. 
We synthesize information from the radio and X-ray maps to assess three wing-formation 
models based on approximate limiting timescales and argue that the hydrodynamic scenario 
faces several difficulties whereas circumstantial evidence favors the reorientation model. Al- 



though 4C +00.58 does not obey the optical-radio correlation of XRGs (ICapetti et al.ll2002 
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Saripalli fc Subrahmanyarul2009l ). confirmation of any of the wing formation scenarios would 



bear on XRG formation generally. 

We thank the referee for helpful suggestions and clarification. E. H.-K. and C S. R. 
thank the support of the Chandra Guest Observer program grants GO90111X and GO89109X. 
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Table 1. Core and Jet Parameters 



Designation 



Distance 

(") 



^1.4 GHz 
(mjy) 



^4.9 GHz 
(mJy) 



-Fx-ray 

(nJy) 



Beq 

(mG) 



Core 


0.0 


41±4 


30±3 


S ci+0.7 


Inner jet 


0.6-6.0 


180±20 


88±9 


1.0 ±0.3 


Inner Knot 1 


2.5 


- 


14±1 


0.3 ±0.2 


Inner Knot 2 


3.3 


- 


15±1 


0.9 ±0.3 


Outer jot 


6.0-10. 


220±20 


70±7 


<0.1 



0.25 


0.8 


0.7±0.1 


- 


- 


0.6 


1.0 


1.1±0.4 


< 7 


7 


0.6=^ 


1.0 


1.0±0.7 


< 8 


9 


0.6=^ 


0.9 


1.5±0.6 


<10 


7 


0.9 


>1.1 


- 


- 


120/5 



^It is not possible to separate the knots at 1.4 GHz, so we use the average value instead of measuring a flux. 

Note. — Most of the jet X-ray emission comes from knot 2 and the outer jet has no X-ray emission. The 
distance is measured radially from the core in arcsec. We report a model X-ray flux from the best-flt power-law 
model with errors reported at 90% confldence. 
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Fig. 1. — : Left: 1.4 GHz VLA A-array map (l."6 x l."3 beam, units in Jy/beam). Right: 
Raw Chandra image (0.3 — 10 keV). The brightest pixel in the AGN has 282 counts. 
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Fig. 2. — : Top Left: 1.4 GHz contours overlaid on smoothed (u = 3 pixels) X-ray image 
with point-like sources excised, clipped at twice the mean background. Top Right: A 
weighted Voronoi tessellation image (S/N = 5.0 in each tile) of the 0.3 — 3 keV events with 
cavities identified. Bottom Left: Negative of an unsharp mask image of top left (using 
a smoothing length of 40 px for subtraction) showing cavities. Bottom Right: Coarsely- 
binned (4x native pixels) image from 0.3 — 5 keV showing extended structures beyond 20 
kpc from the AGN. 
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Fig. 3. — : 1.46 GHz contours overlaid on a combined red+green SDSS image of the host 
galaxy of 4C +00.58. The cyan circle encloses a 20th magnitude resolved "extension" to the 
elliptical galaxy which is slightly bluer. 
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Fig. 4. — : Top: 5 GHz VLA image of the jet with contours (beainsize 0.5 x 0.5", units 
in Jy/beam). Bottom: Smoothed X-ray image with pixel randomization turned off (su- 
perbinned to 1/4-original pixel size) with 5 GHz contours overlaid. 



